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Abstract. Calculations of 7-spcctra for positron annihilation on a selection of 
molecules, including methane and its fluoro-substitutes, ethane, propane, butane 
and benzene are presented. The annihilation 7-spectra characterise the momentum 
distribution of the electron-positron pair at the instant of annihilation. The 
contribution to the 7-spcctra from individual molecular orbitals is obtained from 
electron momentum densities calculated using modern computational quantum 
chemistry density functional theory tools. The calculation, in its simplest form, 
effectively treats the low-energy (thermaliscd, room-temperature) positron as a plane 
wave and gives annihilation 7-spectra that are about 40% broader than experiment, 
although the main chemical trends are reproduced. We show that this effective 
"narrowing" of the experimental spectra is due to the action of the molecular potential 
on the positron, chiefly, due to the positron repulsion from the nuclei. It leads to 
a suppression of the contribution of small positron-nuclear separations where the 
electron momentum is large. To investigate the effect of the nuclear repulsion, as well 
as that of short-range electron-positron and positron-molecule correlations, a linear 
combination of atomic orbital description of the molecular orbitals is employed. It 
facilitates the incorporation of correction factors which can be calculated from atomic 
many-body theory and account for the repulsion and correlations. Their inclusion in 
the calculation gives 7-spcctrum linewidths that are in much better agreement with 
experiment. Furthermore, it is shown that the effective distortion of the electron 
momentum density, when it is observed through positron annihilation 7-spectra, can 
be approximated by a relatively simple scaling factor. 
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1. Introduction 

In the dominant positron annihilation process in matter, a positron annihilates with an 
electron to produce two gamma (7) photons with a total energy of 2mc 2 and a total 
momentum P [1,2]. Nonzero values of P lead to Doppler shifts of the photon energies 
from mc 2 ~ 511 keV. In the case of bound electrons, this Doppler broadening is mostly 
due to the electron momenta [3-5] - the spectra of 7-ray energies (the '7-spectra') 
therefore contain information about the electron momentum distribution in the bound- 
state orbitals and can be used to characterize the orbitals involved [6-9]. 

In light of this fact, positron annihilation spectroscopy is a valuable technique in 
the study of surfaces, interfaces and defects in condensed matter [10,11], of electronic 
properties of quantum dots [12-17], nanoparticles [18,19] and nanocrystals [20], and in 
the study and exploitation of positron annihilation on core electrons [4,8,21-25]. 

In all of these areas, interpretation of the measured annihilation spectra relies 
heavily on theoretical input. However, the problem of calculating 7-spectra for positron 
annihilation on molecules has received little attention to date, with previous positron- 
molecule theoretical studies focussed mainly on understanding the anomalously large 
annihilation rate parameters Z e g [26]. This is in contrast with condensed matter 
systems, where a number of theoretical approaches have been developed [10,11], and 
atomic systems, for which detailed theoretical understanding has been provided by, e.g., 
elaborate variational [5] and many-body theory [6-8] calculations. The aim of this paper 
is to present calculations of 7-spectra for positron annihilation on a range of molecules 
including hydrogen, methane and its fluoro-substitutes, ethane, propane, butane and 
benzene. 

The first measurement of the Doppler-broadened 7-spectrum for positron 
annihilation on molecules was made in 1986 by Brown and Leventhal in low-density 
hydrogen gas [27]. In the early 1990's the confinement of thermalized positrons in a 
Penning trap allowed Tang et al. [28] to measure 7-spectra for a range of molecules, 
including hydrocarbons and perfluorocarbons. This work was built upon by Iwata et 
al. [3] who measured 7-spectra for positron annihilation on a large range of atoms 
and molecules in the gas phase. Despite this extensive set of experimental results, 
there is a paucity of theoretical calculations of the annihilation 7-spectra in molecules. 
Theoretical predictions of the 7-spectra do exist for the molecular ion H^, for which 
(numerical) Coulomb-Born and configuration-interaction Kohn variational calculations 
have been performed [29], and for the small molecules H 2 and N 2 [30,31]. However, 
the application of these sophisticated methods to larger molecules is considerably more 
difficult. More than thirty years ago Chuang and Hogg [32] analysed the annihilation 7- 
spectra for hexane and decane. However, their calculation for the alkane molecules relied 
on analytic wave functions of the carbon orbitals calculated in the 1950's, which are 
relatively crude by modern standards. Furthermore, their calculations disregarded the 
effect of the positron wavefunction. New calculations using modern quantum chemistry 
methods, similar to those reported here, are underway [33]. 
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Although the 7-spectra are determined primarily by the electronic structure, they 
are also sensitive to the positron interaction with the target. Compared to positron 
annihilation on atoms, the positron-molecule problem is significantly more complex 
owing to, e.g., the non-spherical potential experienced by the positron. In this work we 
use two methods to calculate 7-spectra for positron annihilation on molecules. The first 
method is based on electron momentum densities calculated from modern computational 
chemistry tools, such as density functional theory (DFT). This approach was recently 
tested by a number of the authors in a study of positron annihilation on noble gases [9] 
and applied to fluorobenzenes [33]. In its simplest form, this method treats the positron 
as a plane wave and neglects its influence on the 7-spectrum linewidth. To overcome this 
deficiency, and to investigate the effect of the atomic potentials and electron-positron 
correlations, we perform additional calculations that rely on treating the molecular 
orbitals as linear combinations of atomic orbitals. This yields predicted spectra that are 
in good agreement with the computational chemistry calculations and furthermore, it 
gives a framework for the inclusion of the effect of the atomic potentials and correlations 
on the 7-spectrum linewidth (see Ref. [34] for preliminary results). 

The structure of the paper is as follows. Section 2 briefly recaps the general theory 
of positron annihilation 7-spectra. In section 3 the results of the DFT calculations are 
presented for molecular hydrogen, methane and its fluorosubstitutes, ethane, propane, 
butane and benzene. This DFT calculation effectively ignores the effect of the positron- 
molecule interaction on the 7-spectra, and gives 7-spectra that are systematically 
broader than experiment. In section 4 an investigation of the systematic broadening is 
performed. There, we derive the form of the annihilation 7-spectra for molecules using 
a multicentre linear combination of atomic orbital approach. We also show how, in 
this scheme, the effect of the nuclear repulsion and electron-positron correlations can be 
accounted for through momentum-dependent correction factors obtained from positron- 
atom calculations, e.g., using many-body theory. This procedure is applied to H2, CH4 
and CF 4 in section 5. In light of our results, we discuss how the effective distortion of 
the electron momentum density, when it is observed through positron annihilation 7- 
spectra, can be approximated by a relatively simple energy scaling factor. We conclude 
with a summary in section 6. 

2. Theory: gamma-spectra for positron annihilation on molecules 

2.1. Basic equations 

Low-energy positrons annihilate in matter predominately via two-photon production, 
a process in which the total spin of the electron-positron pair must be zero [2,35]. In 
the centre-of-mass frame of the pair, where the total momentum of the two photons 
P = p 7l + p 72 is zero, the two photons propagate in opposite directions and have 
equal energies = p 7 c = mc 2 + \{Ei — Ef) « mc 2 = 511 keV, where Ei and 
Ef are the energies of the initial and final states (excluding rest mass), respectively. 
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When P is nonzero, however, the two photons no longer propagate in exactly opposite 
directions and their energies are Doppler shifted. For example, for the first photon 
E 7l = Ej + mc • V, where V = P/2m is the centre of mass velocity of the electron- 
positron pair and c = cc, where c is the unit vector in the direction of the photon. 
Assuming that ^ < c, and p 7l = E 7l /c ~ mc, the shift of the photon energy from the 
centre of the line, e = E 71 — E 7 , is 

e = mc ■ V = -P • c. (1) 

When a low-energy positron annihilates with a bound electron, the characteristic 
Doppler shifts are determined by the electron energy e n , e ~ ^\e n \mc 2 ^> \e n \. Hence, 
the shift of the line centre E 7 from mc 2 = 511 keV can usually be neglected. The 
annihilation 7-spectrum (Doppler spectrum) can then be written in a form similar to a 
Compton profile, as the integral of the annihilation probability density 

w fl (e,c) = J I%( p )l^( e -^ P - C )^ (2) 

where Afi(P) is the amplitude for annihilation into two photons. It is the fundamental 
quantity to be evaluated. A many-body theory description of it is discussed below. 

For gaseous systems, the measured spectrum represents an average of Wfi(e, c) over 
the direction of emission of the annihilation photons [6, 8] 

1 r°° 

w fi (e) = - p%P)^PdP, (3) 

C J2\e\/c 



where 



is the spherically averaged annihilation momentum density. If the initial state 
corresponds to a positron with momentum k incident on the ground-state target, and 
the possible final states describe the target with an electron missing in orbital n, the 
total spectrum is found by summing over all final states as 

1 f°° 

w k (e) = V u> nk (e) = - / pl{P) AttPoIP, (5) 

where p£(-P) = ^2 n Pnk(P) * s ^e total annihilation momentum density that is equal to 
the sum of the individual annihilation momentum densities of the occupied molecular 
orbitals n. The probability of positron annihilation with core electrons is relatively small 
(owing to the positron repulsion from the nuclei). Hence, the annihilation spectra are 
dominated by the contribution of valence orbitals, except at high Doppler shifts where 
a distinct contribution from the core can be seen [4,8]. 

2.2. The annihilation amplitude 

2.2.1. General form The annihilation amplitude that enters (3) is defined as 

A„ k (P) = (V" 1 J e-^(vMr)£r , (6) 
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where ^(r) and 0(r) are the electron and positron annihilation operators, respectively, 
^k ?+1 i s t ne (fully-correlated) initial state of N electrons and a positron of momentum 
k, and ^l/^ -1 is the (fully-correlated) final state of the target with N — 1 electrons and 
a hole in electron orbital n. It can be calculated using many-body theory methods (see 
Refs. [6-8] for details). 

Figure 1 shows the many-body diagrammatic expansion for the amplitude (see 
the figure caption for an explanation of the diagrams). More specifically, it shows the 
main contributions to the annihilation vertex that describes the short-range electron- 
positron interaction. Figure 1 (a) is the zeroth-order, or independent-particle-model 
(IPM) amplitude, which is given by 

42(P) = / e- tP - r ^(r)^ k (r)rf 3 r, (7) 

where y?k is the wavefunction of the positron, and i/) n is the wavefunction of the 
annihilated electron (hole) in state n. Figures 1 (b) and (c) describe the short-range 
electron-positron correlation corrections to the zeroth-order vertex. Their explicit forms 
can be found in Refs. [6-8]. In addition to the short-range electron-positron correlations, 
there are also correlations between the incident positron and the target, e.g., through a 
polarisation interaction. These interactions are described by a non-local and enegy 
dependent positron self-energy [6-8, 36]. They are included in the fully-correlated 
(or 'dressed') incident positron wavefunction known as the Dyson orbital, which is 
represented by the double line labelled e in figure 1. 




V\ U 2 V\ h>2 V\ v l V\ V'2 




Ml M2 Ml M2 Ml M2 Ml M2 



Figure 1. Main contributions to the annihilation amplitude A„k(P): (a) zeroth-order 
(independent-particle model); (b) first-order correction; (c) r-block (electron-positron 
ladder series) correction. All lines correspond to particles (or holes) propagating 
on top of the ground state of the iV-electron system. The fully-correlated incident 
positron wavefunction (Dyson orbital) is represented by the double line labelled e; the 
wavefunction of the annihilated electron (hole) is labelled by n; the double-dashed lines 
represent the two photons of momentum P; wavy lines represent Coulomb interactions; 
and v and \i are excited positron and electron states, respectively. 



When evaluating the diagrams beyond zeroth order, one must sum over the 
complete set of intermediate electron and positron wavefunctions, including those in 



Effect of positron- atom interactions on the annihilation gamma spectra of molecules 6 



the continuum. This summation can be performed relatively easily for atomic systems, 
since the spherical potential allows for the wavefunctions to be separated into radial and 
angular parts and asymptotic convergence formulae can be used [6-8,36,37]. However, 
the lack of spherical symmetry in the molecular problem makes the development of a 
numerical scheme in which the higher-order diagrams can be calculated a more difficult 
task. | Nevertheless, the shape of the 7-spectra is described to a good first approximation 
by the zeroth-order amplitude [4,6-8], the calculation of which requires only accurate 
positron and electron wavefunctions. We will use this as our starting point. 

To go beyond the zeroth-order approximation and incorporate the effects of the 
higher-order diagrams on the lineshape, we will assume that the positron annihilates in 
the vicinity of one of the atoms of the molecule. The annihilation amplitude diagrams 
can then be calculated using positron and electron orbitals in the field of the atomic 
potentials. 

2.2.2. Plane-wave positron, independent-particle approximation If one assumes that 
the positron is described by a plane wave, then in the low-energy limit k I a.u., one 
has V'k(r) = e lk r ~ 1, for the range of distances at which annihilation occurs. § In this 
case the zeroth-order annihilation momentum density [Eqs. (4) and (7)] reduces to 

PT(P)^ j \MP)\ 2 ^=Pn(n (8) 

where 

MP) = (2n)^' 2 j e" iP ->„(r) d\ (9) 

is the momentum space wavefunction, and p n is the spherically averaged electron 
momentum density of the orbital n.|| The 7-spectrum [cf. (5)] then takes the form 

1 f°° 

w (e) = - p(P) 4wPdP, (10) 

c J2\e\/c 

where p(P) = J2 n Pn{P) is the total electron momentum density of the occupied orbitals 
in the molecule. It is related to the cross section measured using the electron-momentum 
spectroscopy technique [38]. 

2.3. Direct annihilation rate from ^-spectra 

The electron-positron direct annihilation rate in a gas of density n is commonly 
parameterized through the dimensionless quantity Z e g, defined as the ratio of the 
observed annihilation rate A to the rate of free electron-positron annihilation Z c g = 

I In evaluating the diagrams for a molecular system one would have to, e.g., transform the Hartrcc- 
Fock cigenbasis to a (non-orthogonal) one that involves (at least) the atomic orbitals that constitute 
the molecule, centred on their respective atomic sites, as is standard in the derivation of the Roothaan 
equations. Alternatively, one can use a Gaussian basis. 

§ Annihilation occurs at distances of the size of the atom or molecule involved. 
|| The density is normalised by f p n (P) 4irP 2 dP = 1. 
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\jixr\cn [39,40]. By definition, Z e g quantifies the effective number of electrons per 
target atom or molecule, with which the positron can annihilate. In general, it is 
different from the true number of electrons Z owing to the effects of, e.g., the nuclear 
repulsion and positron-electron interactions, which can act to suppress or enhance the 
electron density in the vicinity of the positron. For most molecules, Z e s is also enhanced 
by the vibrational Feshbach resonances [26]. The ratio Z e g/Z therefore gives a measure 
of the strength of such effects. Z eS can be calculated in terms of the Doppler spectrum, 
and hence annihilation momentum density, as 



J-oo J0 

For a plane-wave positron, in the independent-particle-model approximation, it reduces 
to Z c s = Z, the number of electrons in the molecule. 

3. Gamma spectra from the density functional theory calculations 

In this section we present the results of the calculations of the 7-spectra for a selection 
of molecules. They were performed assuming a plane-wave positron, and were based on 
electron-momentum distributions, Eqs. (8)-(10). Specifically, we used the DFT-based 
Becke-3-parameter Lee- Yang-Parr (B3LYP) model [41,42], with the Godbout polarised 
triple-^ valence (TZVP) basis set [43], which is known to produce good agreement 
with the experimental measurements for the molecular orbital momentum profiles [44]. 
These calculations were performed using the computational chemistry package Gaussian 
03 [45]. For C 3 H 8 and C 6 H 6 we used the BP86 functional in GAMESS [46]. The 
calculations employed finite cut-off momentum values in calculating the spectra from 
(10) (in most cases, P max = 3 a.u.). When possible, tails of the form AP~^ were added 
to the orbital momentum densities to obtain more accurate spectra. In most cases this 
did not change their widths by more than 5%. 

Figure 2 shows the results of our plane- wave-positron DFT calculation for the 
annihilation 7-spectra of hydrogen, the fluorosubstituted methanes, and propane and 
benzene. Also shown are the experimental results of Iwata et al. [3] obtained assuming 
that the intrinsic annihilation line shapes can be described by the two-Gaussian form, 



where A, a and b are constants obtained by deconvoluting the intrinsic spectra from the 
measured line shapes. Their values for an extensive range of molecules can be found in 



The values of the calculated and measured full-widths at half-maximum (FWHM) 
are given in table 1. For all of the molecules studied, the calculation systematically 
overestimates the measured line widths, although the main chemical trends are 
reproduced. This can also be seen in figure 3, where we plot the calculated FWHM 
versus the measured linewidths for all of the molecules discussed, and for the noble 




(11) 



10(e) oc exp (-e 2 /a 2 ) + Aexp (-e 2 /b 2 ) , 



(12) 



Ref. [3]. 
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Figure 2. Annihilation 7-spcctra for a selection of molecules: , DFT calculation 

(neglecting the contributions of core orbitals, i.e., C Is and F Is); , DFT spectra 

scaled as described at the end of section 5.1.1; •, experiment (two-Gaussian fit) [3]. 



Table 1. Full widths at half-maximum (keV) of the calculated and measured 7- 
spectra for molecules. 





H 2 


CH4 


CH 3 F 


CH2F2 


CHF 3 


CF 4 


C2H6 


C3H8 


C4H10 




DFT a 


2.14 


2.81 


3.40 


3.81 


4.08 


4.29 


3.05 


2.99 


3.12 


3.14 


Experiment'' 


1.59 


2.06 


2.64 


2.72 


2.71 


2.92 


2.15 


2.18 


2.25 


2.20 



a B3LYP or BP86 calculation of the electron momentum density (i.e., the plane- wave 
positron approximation) without the tails. Core-electron contributions are omitted, 
except for C2H6 and C4H10. 

b Experimental values from the two-Gaussian fit of Ref. [3]. 



gases Ne, Ar, Kr and Xe that were calculated using the Hartree-Fock method [9].^[ 
The graph shows that the FWHM values from plane-wave positron approximation are 
greater than the experimental values by about 40%. 

The main origin of the discrepancy between these calculations and experiment is the 
neglect of the effect of the positron-molecule potential, and in particular, the positron 
repulsion from the nuclei, in the calculation. It leads to an overestimate of the electron- 
positron wavefunction overlap at small positron-nuclear separations, and therefore, of 
the large momentum contributions to the annihilation momentum density, and of the 
large energy shifts in the 7-spectra. 

One should expect that the proper inclusion of the positron-nuclear repulsion 
will lead to a narrowing of the spectra, bringing them more into better agreement 
with experiment. The remainder of the paper outlines a scheme which facilitates this 
inclusion. 

•J Note that for the noble gas atoms, many-body theory calculations are much superior to the 
plane-wave-positron HF method shown here, and produce 7-spectra in excellent agreement with 
experiment [6,8]. 
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Figure 3. FWHM values of the 7-spectra calculated in the plane-wave-positron 
approximation vs. experimental values from Ref. [3], for a range of molecules and 
noble-gas atoms. The dashed line is a linear relation y = 1.4 a;. 



4. Effect of positron-atom interaction 

4-1. Electron momentum density: linear combination of atomic orbitals 

To account for the effect of the nuclear repulsion and electron-positron correlations, it is 
instructive to derive the electron momentum density for a given molecular orbital using 
a multicentre linear combination of atomic orbitals (LCAO). To do this we follow the 
approach of Ref. [47]. 

First, consider a generic molecular orbital ip(r) as a linear combination of the atomic 
orbitals ipi that constitute the molecule, 

^(r) = ^c^(r-Ri), (13) 

i 

where q are the expansion coefficients," 1 " and Rj are the coordinates of the atomic 
centres.* 

Assuming the atomic orbitals can be written in the central-field form ip n im( T ) — 
-Rn«( r )^m(r), where R n i is the radial function and Y im is a spherical harmonic, the 
corresponding momentum space wavefunction (9) is 

Vw(P) = (-i) l R n i(P)Y lm (P), (14) 

+ In general the expansion coefficients are governed by the point-group symmetry of the molecule, and 
by requirements of normalization and energy minimization. 

* In the Born-Oppenheimer approximation the set of the Ri are external parameters as opposed to 
dynamical variables. 
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where R n i(P) = / °° ji(Pr) R n i(r)r 2 dr . The spherically averaged momentum density 
for this atomic orbital, (8), is then given by 

p nl {P) = R 2 nl (P)/4n. (15) 

Hence, the momentum-space wavefunction can then be written as 



A(P) = {-i) k ^pJP)Y lirni {P), (16) 

where i labels the specific atomic orbital in the linear combination (13) and implicitly 
includes the labels for the quantum numbers n and I. The momentum-space molecular 
orbital can then be expressed as follows: 



^(P) = (2tt)- 3/2 J e" iP Xr) d 3 r, (17) 
= ^ Cj e- p - R ^,(P), (18) 

i 

and the corresponding momentum density p(P) is easily shown to be 
p{P) = Y f {-i) ls - U <c jy J Pi {P) Pj {P) j e^^YC n JP)Y hmj (P)dn P . 



(19) 



It is convenient to separate the diagonal and off-diagonal parts in the above double sum, 
and to split the latter into sums involving single-centre and two-centre terms: 

pOp) = En 2 ^( p )+ E pU p )+ E ( 2 °) 

i,j>i 

In this equation, the second term on the right-hand-side is a sum of the products of 
wavefunctions centred on the same atom. They satisfy the usual orthogonality relations 
and take the form 

P%{P) = + cic*)5 hlj 5 mimj ^ Pi (P) Pj (P). (21) 

The third term on the right-hand-side of (20) is a sum over the product of wavefunctions 
centred on different atomic sites. Resolving the plane wave and expressing the integrals 
of three spherical harmonics in terms of 3j symbols (see e.g., [48]), it can be written as 




where R^- = Rj — Rj and we use the notation [I] = (2/ + 1). 

To summarize, in the LCAO approach the molecular orbital momentum density is 
determined by the products of the relevant atomic momentum densities, 



p(P) = Cab, P, Ry)^/pi(P)Pi(P), (23) 
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where C^- are the prefactors given in Eqs. (20)-(22). This quantity represents the 
annihilation momentum density for the corresponding molecular orbital in the low- 
energy plane-wave-positron approximation. A more accurate molecular annihilation 
spectrum will be obtained by replacing the atomic electron momentum densities in (23) 
by the corresponding annihilation momentum densities. 



4-2. Atomic correction factors 

To include the effects of the positron-atom interactions in the molecular annihilation 
momentum density, the plane- wave atomic electron momentum densities pi{P) that 
appear in (23) are multiplied by the atomic correction factor, defined as 



Gi{P) 



(24) 



Pi(py 

i.e., as the ratio of the true annihilation momentum density of the atomic orbital i to 
its plane-wave independent-particle approximation. 

For atomic hydrogen, the correction factors can be calculated accurately from many- 
body theory (see figure 4 (a) for a diagrammatic definition). They therefore include the 
effect of electron-positron and positron-atom correlations, in addition to that of the 
nuclear repulsion. A simpler correction factor that can be used for any atom and still 
include the effect of the atomic potential and nuclear repulsion can be obtained using 
a Hartree-Fock (HF) positron function with the zeroth-order (independent-particle) 
annihilation vertex (see figure 4(b)). Note that the effect of the short-range electron- 
positron correlations are not included in this correction factor. 
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(a) (b) 

Figure 4. Atomic correction factors: (a) the 'exact' atomic correction factor for 
atomic orbital n, defined as G(P) = p„(P) I Pn(P) , where is calculated using the 
positron Dyson orbital and the full vertex (see figure 1); (b) the simpler Hartree-Fock 
correction factor defined as G^ F (P) = p^' HF (P)/p n (P), where p^' HF (P) is calculated 
with a Hartree-Fock positron wavefunction and the zeroth-order annihilation vertex. 



We will apply this procedure to calculate the 7-spectra of H 2 , CH 4 and CF 4 . Figures 
5 and 6 show the atomic correction factors for hydrogen, carbon and fluorine, calculated 
for the s-wave positron wavefunction which dominates at the low (room-temperature) 
momentum k = 0.04 a.u. To highlight the effect that these factors have on the shape of 
7-spectra, they are shown normalised to unity at P = 0. 
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Momentum P (a.u.) 

Figure 5. Momentum-dependent atomic correction factors for hydrogen calculated 
with the Hartree-Fock (i.e., static field) and Dyson positron orbitals using different 
approximations for the annihilation vertex. All correction factors are normalised to 
unity at P = 0. The oscillations in the tails are artefacts arising from computational 
"noise" in the many-body theory calculations. 

Figure 5 shows the correction factor for hydrogen calculated using different 
approximations to the annihilation vertex and positron wavefunction. The effect of 
the nuclear repulsion is most easily seen through the zeroth-order factor G^ F . The 
shape of this factor is mainly governed by the nuclear repulsion. It is evident that 
the effect of the nuclear repulsion on the positron is to suppress the high-momentum 
regions of the momentum density. Further suppression results from inclusion of the 
additional correlations. One can also see that for a given approximation to the positron 
wavefunction, adding the first-order and higher-order (r-block) correlation corrections 
to the vertex increases the degree of suppression of high momentum components to 
the annihilation momentum density. Physically, this is explained by noting that the 
corrections to the vertex involve positron annihilation with spatially diffuse virtual 
electrons. On the other hand, using the correlated positron Dyson orbital instead of 
the Hartree-Fock (i.e., static field) positron wavefunction increases the high- momentum 
component for a given approximation for the vertex. Again this is what one should 
expect physically, since the inclusion of the attractive correlational potential accelerates 
the positron towards the atom and increases the momentum of the annihilating electron- 
positron pair. Regarding the 7-spectrum shape therefore, there is some compensation 
between including the correlation corrections to the vertex and increasing the accuracy 
of the positron wavefunction. However, the dominant effect of G(P) is always the 
suppression of high momenta. 

For the carbon and fluorine atoms, we have only the zeroth-order Hartree-Fock 
correction factors at our disposal}! and these are shown for the ground states orbitals in 

jj The many-body theory cannot be applied directly to open-shell atoms, but the static-field effect is 
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Figure 6. Momentum-dependent atomic correction factors for the ground state 
orbitals of hydrogen, carbon and fluorine calculated using the Hartree-Fock positron 
wavefunction with the zeroth-order annihilation vertex. All correction factors are 
normalised to unity at P = 0. 



figure 6. Note that the correction factors for the 2s orbitals show singularities arising 
from the nodes in the Hartree-Fock wavefunctions (which result in the zeros of the 
corresponding electron momentum densities p2 S {P))- Nevertheless, the suppression 
of high momentum contributions to the spectra is clear. Moreover, when evaluating 
the 7-spectra in the LCAO approach we circumvent the problem of the singularities 
in the correction factors by simply replacing the atomic electron momentum densities 
with the annihilation momentum densities, i.e., we calculate p"(-P) directly rather than 
Pi(P)Gi(P). For carbon and fluorine, the positrons predominatly annihilate on the 
valence electrons. The magnitude of the correction factor for the valence electrons is 
therefore substantially greater than for the core electrons. 

Finally, although for hydrogen the correction factors are sensitive to the positron- 
atom and positron-electron correlations, for heavier atoms most of the suppression at 
high momenta is due to the effect of the nuclear repulsion. Hence, for such atoms the 
HF correction factor will describe the true suppression more accurately than it does for 
hydrogen. This point will be discussed further in setion 5. 

5. Gamma spectra with atomic correction factors 

5.1. Molecular hydrogen 

5.1.1. Gamma spectrum of H2 In the LCAO approximation, the ground state 
molecular orbital ij) a {v) can be written a linear combination of two hydrogenic Is 

the dominant feature for higher- atoms and core orbitals. 
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Momentum P(a.u.) 

Figure 7. Hydrogen molecule electron momentum densities calculated using the 
B3LYP/TZVP DFT (black solid curve) and LCAO, (26) (dots), and annihilation 
momentum densities obtained from (26) by multiplying the "plane-wave" (PW) 
atomic density pi s (P) by the correction factors G(P) obtained from many-body-theory 
calculations using various approximations to the annihilation vertex and positron 
wavefunction (see legend). All densities are normalised to unity at P = 0. 



orbitals, each centred on one of the protons, 

^ 9 (r) = — 1 hfc.(r - R x ) + <Mr - R 2 )] , (25) 

V A 1 + b) 

where ^>i s (r) = {jice')~ l ^ 2 e~ ar and S is the overlap integral f ipi s (r — R 1 )^ ls (r — 
R 2 )d 3 r. It can be calculated analytically using prolate spheroidal coordinates as 
S = e~ aR [l + aR + |(ai?) 2 ] , where R is the internuclear distance [49]. A variational 
calculation gives a ~ 1.19 for R = 1.38 a.u., and therefore, S ~ 0.69 [50]. 

Substituting (25) into (20) gives the molecular orbital momentum density 

= T^s pls(P) V- +j °( PR ^ ' ( 26 ) 

where pi s (P) is the momentum density of the hydrogenic Is orbital, 

»u{P) = ^£pr (27) 

Figure 7 shows that the LCAO molecular orbital density given by (26) is very close to 
that from B3LYP /TZVP density functional theory calculation. 

Figure 7 also shows the hydrogen molecule annihilation momentum densities 
obtained from (26) by scaling the atomic momentum density with the correction factors 
G(P) obtained from many-body-theory calculations using various approximations to 
the annihilation vertex and positron wavefunction (cf. figure 5), i.e., by substitution 
Pn(P) — > G(P)pis(P). To emphasise the effect of the correction factors on the 
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Energy shift (keV) 

Figure 8. Annihilation 7-spectra of H2 calculated using LCAO in the plane-wave- 
positron approximation, i.e., using the density from (26) (long-dashed line), and 
incorporating the 'exact' momentum-dependent atomic correction factor G Z\, r (P) 
(solid line) and the energy-shift-dependent correction factor G 7 (e) (crosses, see text); 
solid circles show the measured spectrum, as represented by the two-Gaussian fit, 
Ref. [3]. All spectra are normalised to unity at e = 0. 



annihilation momentum density, all densities are normalized to unity at P = 0. As 
expected, the use of correction factors leads to significant narrowing of the momentum 
ditribution. 

The corresponding annihilation 7-spectra are shown in figure 8, where they are 
compared with the experimental result of Iwata et al. [3]. The FWHM values of the 7- 
spectra are given in table 2. The spectra obtained from the B3LYP/TZVP and LCAO 
electron momentum densities are practically indistinguishable, and we show only the 
latter in the figure. This 'low-energy plane-wave positron' calculation gives the FWHM 
value about 35% greater than experiment. As discussed above, this overest filiation of 
the high-momentum contribution is largely due to the neglect of the positron-nuclear 
repulsion. Introducing the effects of the nuclear repulsion through the zeroth-order 
vertex Hartree-Fock correction factor G^ F , narrows the spectrum, but it remains ~ 18% 
broader than experiment. Using the most accurate, 'exact' atomic adjustment factor 
G^J_ B 1+T (P) gives the FWHM that is in very good agreement (within 6%) of the measured 
value. This shows that for H 2 , the inclusion of the nuclear repulsion provides only 
50% of the required narrowing to bring the plane-wave calculation into agreement with 
experiment, with the remaining 50% being due to the electron-positron correlations. 
However, as the nuclear potential becomes stronger in molecules containing heavier 
atoms, the effect of the potential on the linewidth is expected to become the principal 
cause of the narrowing, overriding the effect of the correlations. 
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Table 2. FWHM of the annihilation 7-spectra and annihilation rate parameter Z e ff 
for H2 calculated in the plane- wave positron approximation (B3LYP/TZVP DFT and 
LCAO) and using the atomic correction factors in different approximations, and results 
of other theoretical calculations and experiments. 



Method 


FWHM 

(keV) 




B3LYP/TZ-DFT (plane- wave) a 


2.14 


2.00 


LCAO (plane-wave) a 


2.14 


2.00 


LCAOxG^ F (P) a 


1.88 


0.78 


LCAOxG^f^F) 3 - 


1.73 


1.98 


LCAOxG? F , , JP) a 


1.40 


4.59 


LCAOxG^ ys (P) a 

1 J V 1 IV / s\ v * r\ \-Lj 


2.01 


2.62 


LCAOxGo|i(P) a 


1.85 


6.54 


LCAOxG°f; +r (P) a 


1.50 


13.9 


Kohn variational 13 




10.3 


Kohn variational 




9.8 


Kohn variational 01 




12.6 


Stochastic variational 




15.7 


Molecular iJ-matrix 




10.4 


Expt. g 


1.59 




Expt.' 1 




14.7 


Expt. 1 




16.0 


Expt. j 




14.6 


a This study, R = 1.4 a.u., k = 0.04 a.u. 
b 'Method of models' calculation from Ref. [51]. 
Calculation for R=1A a.u., k = 0.04 a.u., Ref. [5 
d 'Method of models' calculation from Ref. [53]. 



Confined variational method calculation, Ref. [54]. 
f Molecular R- matrix with pseudo-states, Ref. [55]. 
s Experiment, room-tcmpcrature positrons, Ref. [3]. 
h Experiment, room-temperature positrons, Ref. [56]. 

1 Experiment, room-temperature positrons, Ref. [57]. 
1 Experiment, room-temperature positrons, Ref. [58]. 

In addition, we can test a simpler method in which the correction factor is applied 
directly to the Doppler-shift spectrum rather than in the two-7 momentum density. The 
method is suggested by an observation that the momentum densities and the 7-spectra in 
figures 7 and 8 are approximately Gaussian (i.e., parabolic in the semilogarithmic plots). 
If one assumes that both the electron momentum and the true annihilation momentum 
atomic densities are Gaussians, then then the atomic correction factor G(P), (24), will 
also be a Gaussian. It is easy to see from Eqs. (5) and (10) that in this case the 
corresponding annihilation spectra are also Gaussian. Their widths are proportional to 
those of the momentum densities (scaled by c/2). Denoting the FWHM values of the 
true annihilation spectrum and that found using the electron momentum density by a a 
and a m , we can introduce the atomic 7-spectrum correction factor 6? 7 (e) = w(e)/w m (e), 
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which will also be a Gaussian, of width a G = (cr~ 2 — cr~ 2 ) 1 ^ 2 . 

This means that if the 7-spectrum of molecular hydrogen calculated using the 
electron momentum density (i.e., plane- wave positron) is close to a Gaussian, of width 
<7^ 2 , then one can approximate the true annihilation 7-spectrum of H2 by u> H2 (e) ~ 
G 1 {e)w i ^{e). The latter will also be a Gaussian with the width 

-1/2 a H 2 

(28) 



a? 2 



{el 2 ) 2 °l 



o- 



m 2 /o- G y 



Figure 8 shows that adjusting the the spectrum by the Doppler-shift-dependent 
correction factor G 1 (e) for H2 produces a spectrum in good agreement with the result 
obtained using G(P) and with experiment. 

Estimating the width of G y (e) from the best calculation shown in figure 6 (G^.p) 
gives oq = 2.10 keV. Using this value together with the plane- wave-positron width 
<t^ 2 = 2.14 keV (table 2) in (28) gives the annihilation spectrum width crj^ 2 ~ 1.50 keV, 
in agreement with the best calculated value in table 2. 

This example shows that oq a m , i.e., that the width of the positron atomic 
correction factor is close to the that of the electron momentum distribution. In this case 
the expression under the square root in (28) is close to 2, giving a a ~ a m /lA. It appears 
that this relation is quite general, as shown by the FWHM values for other molecules 
and atoms in table 1 and figure 3. This suggests that to obtain a closer description of 
experiment, the spectra calculated neglecting the positron can be modified by simply 
scaling the energy shift by a factor 1.4. The corresponding spectra are shown in figure 2 
by dotted lines. The agreement with experiment is excellent for hydrogen and good for 
alkanes (except in the wings), with larger discrepancies for molecules containing fluorine 
atoms. 



5.1.2. Z c g of H 2 The correction factor approach also provides a means of estimating the 
normalized annihilation rate parameter Z e g- for positrons on H 2 . Its accurate theoretical 
calculation has proved to be a difficult problem to which considerable attention has 
been paid. Extensive calculations have been performed by Armour and co-workers, 
who have applied various versions of the Kohn variational method [51-53]. However, 
their best theoretical results still underestimate commonly accepted measured value of 
14.6 [58]. Only recently has this problem been effectively solved, with Zhang et al. [54] 
performing a stochastic variational method calculation which provided an accurate 
total wavefunction of this two-elelctron-one-positron system and yielded a low-energy 
(thermalized) Z e g of 15.7. 

In this work we calculate the Z e g from the 7-spectrum using (11). The results of 
the B3LYP/TZVP DFT and LCAO calculations using the atomic correction factors are 
given in table 2. As one should expect, the plane- wave positron approximation (both 
B3LYP/TZVP and LCAO) gives Z e g = 2, i.e., the number of electrons. The effect 
of the repulsive atomic field, implicit in G^ F , reduces the magnitude of Z e g by more 
than a factor of two relative to the plane-wave calculation, owing to the reduced overlap 
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between the electrons and the incident positron. However, the successive inclusion 
of higher-order correlation effects through the correction factors produces a marked 
increase in the calculated Z e ff, with the full many-body theory (Dyson positron orbital, 
exact vertex) giving Z e g = 13.9, which is within 5% of the measured value. This level 
of agreement is partly fortuitous, but it is still remarkable that the calculation which 
uses atomic correction factors gives a molecular Z c q value in such close agreement with 
experiment. The fact that it does gives some a posteriori justification for the LCAO 
correction factor approach and its ability to take appropriate account of the positron- 
molecule and positron-electron correlations in H 2 . 



5.2. Gamma spectrum of methane 

Methane, CH4, has a tetrahedral structure with a point group The four hydrogenic Is 
orbitals span the reducible representation r^. ed = Ai®T 2 , where Ai is the fully symmetric 
one-dimensional representation and T2 is a degenerate three-dimensional representation. 
Denoting the carbon ground state orbitals by ip n i and the four atomic hydrogen ground- 
state orbitals by <7j, the ground state molecular LCAO of methane are (see e.g., [59]) 

if> lai (r) = a^u + M 2s + 7l a Al , (29) 
^ 2ai (r) = a 2 t/j ls + f3 2 ip 2s + 72 a Al , (30) 

V'lta.,,,, 0) = ^2 Px , y ,z + (31) 

where a, (3, 7 and 5 are the expansion coefficients and the symmetry adapted linear 
combinations (SALC) of pendant atoms are given by 

a Al (r) = o\ + a 2 + <t 3 + a 4 , 

aj 2 (r) = cti - o 2 - 03 + a 4 , 

a T y 2 (r) = cti + ct 2 - 0-3 - 0"4, 

o-J 2 (r) = 01 - cr 2 + o- 3 - 4 - 

The total electron momentum density for methane is given by the sum of the 
respective individual molecular orbital densities, 

p(P) = 2p lai (P) + 2p 2ai (P) + 6 Plt2 (P) , (32) 

obtained from Eqs. (29)-(31) and (20) in terms of the atomic orbital momentum 
densities, as 



Pnai (P) = a 2 nP c ls {P) + p 2 nP g(P) + 2a n f3 n ^ P f s (P)pg(P) + 4^(P) [1 + 3 j (PP HH )] 
+8 Jo (PR CH ) ln J P l(P) \aJpUP) + Pn\/pg(P)\ , (33) 



Plt2 (P) = 5 2 P2 P (P) + HpUP) [1 " Jo(PRur)] + &y a fij 1 {PRo B )y/f$,{P)ft,{P). (34) 

where Pch an d Phh are the C-H and H-H bond lengths, respectively. 

We construct the molecular orbitals of Eqs. (29)-(31) in two ways. In the first 
case, we use Slater-type atomic orbitals (LCAO-STO) optimized for the molecular 
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Figure 9. Electron momentum densities (i.e., plane- wave positron annihilation 
momentum densities) for the molecular orbitals of methane: (solid lines) 
B3LYP/TZVP; (dashed lines) Slater-type orbitals; (dotted lines) Hartree-Fock. 

environment [60]. In the second case, we use Hartree-Fock atomic orbitals (LCAO- 
HF). In both cases the LCAO expansion coefficients are taken to be those of Ref. [60], 
with a C-H distance of Rqh = 2.05 a.u. 

Figure 9 shows the electron momentum densities calculated using the DFT based 
B3LYP/TZVP model, and the LCAO result from (32) using the Slater-type and Hartree- 
Fock orbitals. For a given molecular orbital, the results of all three calculations are 
similar. However, in general the LCAO-STO calculation is in better agreement with 
the B3LYP/TZVP calculation than the LCAO-HF one, owing to the fact that the HF 
orbitals used are purely atomic in nature, i.e., they are in no way adjusted to account 
for the molecular environment. ff Note that the momentum density of the lai orbital 
(predominantly carbon Is in nature) is much broader than those of the valence 2cti and 
lt 2 orbitals due to the fact that the core electrons have larger characteristic momenta. 

To take into account the effect of the positron-atom interactions on the annihilation 
momentum density and 7-spectrum, we replace the atomic electron momentum densities 
in Eqs. (33) and (34) by the corresponding annihilation momentum densities calculated 
using the Hartree-Fock positron wavefunction. These are then used in Eq. (32) and the 
annihilation spectrum is found from Eq. (5). This procedure is equivalent to the use of 
atomic correction factors G^ F (P), which is expected to narrow the spectrum (see figure 
6). 

Figure 10 shows the annihilation 7-spectra of CH 4 , calculated using the 
B3LYP/TZVP and LCAO (both STO and HF), and the LCAO-HF results corrected for 
the positron-atom interactions, as well as the measured spectrum. The B3LYP/TZVP 

ft This is in contrast with the treatment of H2 in section 5.1, where we used hydrogenic orbitals with 
a = 1.19 from a variational calculation, rather than those of the H atom (a = 1). 
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Figure 10. Annihilation 7-spectra of CH4: (dashed line) B3LYP calculation; (dotted 
line) LCAO-STO; (dot-dashed line) LCAO-HF (all correspond to the low-energy plane- 
wave positron); (thick solid) LCAO-HF using the atomic annihilation momentum 
densities; (circles) experiment, two-Gaussian fit, Rcf. [3]. All spectra are normalised 
to unity at e = 0. 

spectrum is significantly broader than the measured one, as seen earlier in figure 2. The 
fact that the LCAO-STO result is in very good agreement with the B3LYP/TZVP result 
confirms the applicability of the LCAO approach. The LCAO-HF spectrum calculated 
in the same low-energy plane- wave positron approximation (i.e., neglecting the positron 
wavefunction) is narrower than B3LYP/TZVP, but still much broader than experiment. 
The overestimation of the width must therefore originate from the neglect of the positron 
wavefunction. Note that the plane-wave calculation overestimates the spectral FWHM 
significantly, even when the la\ core orbital is neglected (see table 3). 

Correcting the LCAO-HF spectrum by using atomic annihilation momentum 
densities in Eqs. (32)-(34) produces a 7-spectrum which is very close to experiment 
for Doppler shifts smaller than 3 keV (solid line in figure 10). The excellent agreement 
with experiment is somewhat fortuitous because the starting point for the adjustment, 
i.e., the plane- wave LCAO-HF calculation, does not perfectly agree with the accurate 
B3LYP/TZVP molecular calculation. However, the need for and effect of the correction 
is clear. As can be seen in table 3, a considerable narrowing is observed with the FWHM 
reduced from 2.52 keV to 2.09 keV, compared with the measured value of 2.06 keV. 
Thus most of the required narrowing to obtain agreement with experiment has been 
accounted for through the effective inclusion of the nuclear repulsion with all of the 
atoms. In the case of H 2 we saw that electron-positron annihilation vertex corrections 
further narrow he annihilation spectrum (see figure 5). Therefore, it is expected that if 
one combined more accurate atomic correction factors with the broader LCAO-STO (or 
B3LYP/TZVP) electron momentum density (see figure 10), the resulting annihilation 
spectrum would provide a good description of the experimental data. 
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Table 3. FWHM of the 7-spectra for CH 4 (in keV). 



Orbital 


B3LYP 




LCAO 




Exp. 






STO 


HF 


HF a 






1.96 


1.85 


1.64 


1.45 




1*2 


3.31 


3.52 


2.98 


2.55 




Total (valence) 


2.81 


2.85 


2.43 


2.09 




Total (all) 


2.97 


2.97 


2.52 


2.09 


2.06 



a Corrected by using the atomic annihilation momentum densities. 

Finally, note that a broad core contribution can be seen clearly in the experimental 
spectrum for Doppler shifts > 3 keV. To reproduce this feature, the theory would need 
to be able to calculate accurately the size of the core contribution relative to that of the 
valence molecular orbitals [8]. 

5.3. Gamma spectrum of tetrafluoromethane 

Like methane, tetrafluoromethane (CF 4 ) belongs to the point group T^. Similar to the 
hydrogens in CH4, the groups of fluorine Is, 2s and 2p z orbitals (with the z axis directed 
locally towards the carbon atom) span the reducible representation T™ d = A\ © T 2 . 
In contrast to CH 4 , for CF 4 one has the additional complication of tt bonding, with 
the complete set of the eight fluorine 2p x y orbitals forming a basis for a reducible 
representation r; cd = E®T 1 @T 2 (see e.g., [50,59]). Accordingly, the B3LYP/TZVP 
calculation for CF 4 gives the following electronic orbital structure: core + laf lt\ 2a\ 
2t\ e 4 3^2 tf, where the core consists of a triply degenerate T2 symmetry orbital and 
a fully symmetric orbital that are predominantly F Is in nature, and another fully 
symmetric orbital which is predominantly C Is in nature. Owing to the repulsion from 
the nuclei, the positron will have difficulty reaching and subsequently annihilating on 
these core orbitals, and the effect of the correction factors will be to greatly suppress 
the magnitude of the contribution of these core orbitals in comparison with the valence 
orbitals. 

To employ the momentum- dependent correction factors, note that the e and t\ 
molecular orbitals are nonbonding, and are composed entirely of F 2p orbitals. Their 
respective electron-momentum densities can therefore simply be multiplied by Gp2p{P) 
when calculating the corrected annihilation spectra. To apply the correction to the 
remaining molecular orbitals, however, we must construct the LCAO molecular orbitals 
and replace the electron momentum densities with the annihilation momentum densities 
calculated, in this case, in the Hartree-Fock zeroth-order vertex approximation, as was 
done for CH 4 . The molecular orbitals are 

l/j ai = atl/jcis + Mc2s + 7l0p£i + 5 1°F2 S + Vl°F2p x 

where a, fl, 7, S and rj are the expansion coefficients, and the a and ir SALC orbitals 
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Figure 11. Electron momentum density of the valence orbitals of CF4 calculated 
using the B3LYP (solid lines) and the LCAO-HF approach (dashed lines). Dot-dashed 
lines show the annihilation momentum densities obtained in the LCAO-HF approach 
using atomic orbital correction factors. 



are given in Ref. [61]. 

Figure 11 shows the momentum densities for the valence molecular orbitals 
calculated in the LCAO approach from (20) using the Hartree-Fock atomic orbitals. 
We determined the coefficients in the LCAO expansions by fitting the LCAO-HF 
momentum densities to the accurate B3LYP results (see figure 11), using the coefficients 
from Ref. [61] as the starting approximation. This approach yields the LCAO-HF 
electron momentum densities in good agreement with the B3LYP results, especially 
for the t2 orbitals which contain most of the electrons. The total valence electron 
momentum density obtained in the LCAO-HF approach is very close to that of the 
B3LYP calculation. Also shown in figure 11 are the annihilation momentum densities 
obtained in the LCAO-HF approach using atomic annihilation momentum densities 
from the zeroth-order Hartree-Fock calculation (i.e., including the G^ F (P) factors for 
the atomic orbitals involved). As expected, the correction factors reduce the magnitudes 
of the momentum densities (because of the positron repulsion from the nuclei) and, in 
particular, suppress the contribution of high momenta. 

Figure 12 shows the calculated annihilation 7-spectra along with the the measured 
spectrum. The corresponding FWHM values are given in table 4. The spectra from the 
B3LYP and LCAO-HF calculations, which neglect the positron-molecule interaction, are 
in good agreement with each other, but are significantly broader than the experiment. 
On the other hand, the LCAO-HF calculation that incorporates the atomic correction 
factors gives results much closer to the measured spectrum. For CF 4 , the dominant 
contribution to the annihilation momentum density comes from the twenty fluorine 2p 
electrons. Therefore, in addition, the figure shows the result of a simplified LCAO-HF 
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Figure 12. Annihilation 7-spectra of CF4: (long-dashed line) B3LYP calculation; 
(short-dashed line) LCAO-HF (both correspond to the low-energy plane- wave positron 
approximation); (solid line) LCAO-HF with atomic correction factors; (crosses) LCAO- 
HF with the fluorine 2p-orbital correction factor; (circles) experiment, two-Gaussian 
fit, Ref. [3]. The spectra were calculated from the total densities shown in figure 11 
with power-type momentum density tails, and normalised to unity at e = 0. 

Table 4. FWHM of the 7-spectra of CF 4 (in keV). 



B3LYP 




LCAO 




Exp. 




HF 


HF a 


HF b 




4.48 


4.52 


3.32 


3.17 


2.92 



a Corrected by using the atomic annihilation momentum densities. 
b Corrected by using the fluorine 2p-orbital factor Gf2 P (-P)- 

calculation in which we use the correction factor of the fluorine 2p atomic orbital for 
all molecular orbitals. The results of both adjustment methods are very similar. This 
suggest that if a specific atomic orbital dominates the molecular annihilation spectra, 
then the spectra can be brought in line with experiment using only the correction factor 
of this orbital. 

The inclusion of the correction factors causes ~ 30% reduction in the linewidth, as 
documented in table 4, and brings it into improved agreement with the measured value. 
Note that one can gauge the relative importance of the inclusion of the atomic potential 
and the inclusion of the proper vertex in fluorine by considering these effects in neon. For 
the valence 2s and 2p orbitals in neon, the plane- wave positron linewidth of 4.94 keV [9]. 
This compares with the FWHM of 3.82 keV for the HF zeroth-order-vertex calculation, 
and 3.55 keV, when using the Dyson orbital and full 0th + lst-order +T vertex [6,62]. 
Thus for neon the majority of the narrowing comes from the inclusion the effect of the 
atomic potential, and we should also expect this to be the case for CF 4 . The zeroth- 
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order HF correction factor therefore brings the calculation closer to the measured value 
than was observed for the corresponding quantity in H 2 , where the nuclear potential is 
not as strong. 

6. Summary and conclusions 

Two approaches have been used to calculate 7-spectra for positron annihilation for a 
selection of molecules including methane and its ffuorosubstitutes. The first method 
used modern quantum chemistry tools and density functional theory (B3LYP or BP86) 
to calculate the electron momentum densities, which were then converted into the 
annihilation spectra. This method assumes that a positron thermalized at room 
temperature can be described using the plane- wave approximation, i.e., neglecting its 
interaction with the target. The corresponding annihilation spectra were found to be 
systematically broader than experiment by about 40%. 

These calculations were complemented in an approach that relied on expanding 
the molecular orbitals as linear combination of atomic orbitals. It allowed us to 
incorporate the effect of positron-atom interactions and electron-positron correlations 
using momentum-dependent atomic correction factors. In particular, this method 
enables one to include the effect of positron-nuclear repulsion, which leads to a significant 
narrowing of the annihilaton spectra compared to the plane-wave approximation. 

For the hydrogen atoms the atomic correction factor could be calculated near 
exactly from many-body theory. Its inclusion produced a 7-spectrum of H 2 in good 
agreement with experiment. Furthermore, for H2 the method gave the absolute 
annihilation rate parameter Z c g rj 14 for thermal room-temperature positrons. This 
value is close to the commonly accepted experimental Z e g = 14.6 [58] and best 
variational calculation, Z eS = 15.7 [54]. 

For carbon and fluorine we used much simpler, Hartree-Fock positron correction 
factors. Nevertheless, the effective inclusion of the nuclear repulsion through these 
correction factors resulted in a significant narrowing of the spectra of CH 4 and CF 4 , 
giving much improved agreement with experiment. For a polyatomic molecule the 
annihilation spectrum is determined by contributions of many molecular orbitals. In 
CF 4 the dominant orbitals are those that are fluorine- 2p in nature, which allows one to 
adjust the spectrum using a single atomic correction factor. The remaining discrepancies 
are due to the neglect of the positron-atom and positron-electron correlations. 

Any complete ab initio theory of positron annihilation on molecules must include 
the full electron and positron dynamics, accounting for the positron-nuclear repulsion 
and electron-positron correlation effects. However, we have shown that realistic 
annihilation 7-spectra for molecules can be obtained from a standard LCAO electron 
momentum density approach, by replacing the atomic electron momentum densities 
(equivalent to plane-wave annihilation momentum densities) with atomic annihilation 
momentum densities. 

We have shown that the effective distortion of the electron momentum density, 
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when it is observed through positron annihilation 7-spectra, can be approximated by a 
simple scaling which reduces the momenta and 7-ray Doppler shifts by a factor of 1.4. 
The origin of this scaling is in the similarity between the electron orbital momentum 
densities and corresponding momentum-dependent positron correction factors. Realistic 
annihilation 7-spectra can therefore be obtained from standard electron momentum 
density calculations by including this scaling. Conversely, the simple scaling factor 
allows for the deduction of electron momentum densities from positron annihilation 
7-spectra. 
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